A Fast Cyclic Reduction Algorithm for A Quadratic Matrix Equation Arising from Overdamped Systems
نویسندگان
چکیده
We are concerned with a class of quadratic matrix equations arising from the overdamped mass-spring system. By exploring the structure of coefficient matrices, we propose a fast cyclic reduction algorithm to calculate the extreme solutions of the equation. Numerical experiments show that the proposed algorithm outperforms the original cyclic reduction and the structure-preserving doubling algorithm. Keywords—Fast algorithm, Cyclic reduction, Overdamped quadratic matrix equation, Structure-preserving doubling algorithm
منابع مشابه
Detecting and Solving Hyperbolic Quadratic Eigenvalue Problems
Hyperbolic quadratic matrix polynomials Q(λ) = λ2A + λB + C are an important class of Hermitian matrix polynomials with real eigenvalues, among which the overdamped quadratics are those with nonpositive eigenvalues. Neither the definition of overdamped nor any of the standard characterizations provides an efficient way to test if a given Q has this property. We show that a quadratically converg...
متن کاملOn the solution of algebraic Riccati equations arising in fluid queues
New algorithms for solving algebraic Riccati equations (ARE) which arise in fluid queues models are introduced. They are based on reducing the ARE to a unilateral quadratic matrix equation of the kind AX2 + BX + C = 0 and on applying the Cayley transform in order to arrive at a suitable spectral splitting of the associated matrix polynomial. A shifting technique for removing unwanted eigenvalue...
متن کاملAlgorithms for hyperbolic quadratic eigenvalue problems
We consider the quadratic eigenvalue problem (QEP) (λ2A+λB+ C)x = 0, where A,B, and C are Hermitian with A positive definite. The QEP is called hyperbolic if (x∗Bx)2 > 4(x∗Ax)(x∗Cx) for all nonzero x ∈ Cn. We show that a relatively efficient test for hyperbolicity can be obtained by computing the eigenvalues of the QEP. A hyperbolic QEP is overdamped if B is positive definite and C is positive ...
متن کاملOn iterative methods for the quadratic matrix equation with M-matrix
In this paper, we consider Newton’s method and Bernoulli’s method for a quadratic matrix equation arising from an overdamped vibrating system. By introducing M-matrix to this equation, we provide a sufficient condition for the existence of the primary solution. Moreover, we show that Newton’s method and Bernoulli’s method with an initial zero matrix converge to the primary solvent under the pro...
متن کاملA New Two-stage Iterative Method for Linear Systems and Its Application in Solving Poisson's Equation
In the current study we investigate the two-stage iterative method for solving linear systems. Our new results shows which splitting generates convergence fast in iterative methods. Finally, we solve the Poisson-Block tridiagonal matrix from Poisson's equation which arises in mechanical engineering and theoretical physics. Numerical computations are presented based on a particular linear system...
متن کامل